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^ '. Abstract 

p • We analyze the branching ratio for the FCNC mode -n^uij in the stan- 

^ \ dard model with QCD effects taken into account consistently to next-to-leading 

order. This involves a two-loop renormalization group analysis for the charm 
^ ■ contribution, presented in this paper, and the calculation of 0{as) corrections to 

\ all orders in mt/M^ for the top-quark case that we have described elsewhere. 

The inclusion of next-to-leading corrections reduces considerably the theoretical 
uncertainty due to the choice of the renormalization scales, inherent in any calcu- 
lation to finite order in perturbation theory. For n^uu this point has not 
been discussed previously. In particular, the related uncertainty in the determi- 
nation of \Vtd\ from B{K^ tt^uu) is reduced from ~ 30% to ~ 7% for typical 
values of the parameters. Simultaneously also the dependence of B{K~^ —>■ n^uu) 
on the choice of mc is considerably reduced. We also give the next-to-leading order 
expression for the short- distance part of Kl — > Impacts of our calculations 
on the determination of the unitarity triangle are presented. 



* Supported by the German Bundesministerium fiir Forschung und Technologic under con- 
tract 06 TM 732 and by the CEC science project SC1-CT91-0729. 



1 Introduction 



The rare decay mode n^uu is one of the most interesting and promising 

processes in the rich field of rare decay phenomenology ^ and has already 
received considerable attention in the past |, H, |^ . The main reasons for 
this interest may be summarized as follows: 

— i> n^uu is a "clean" process from the theoretical point of view: 

• It is short- distance (SD) dominated. Possible long-distance contributions 
have been considered in |TD[ and found to be negligibly small (three orders 
of magnitude smaller than the SD contribution at the level of branching 
ratios) . 

• Being a semi-leptonic process, the relevant hadronic operator is just a cur- 
rent operator whose matrix element can be extracted from the leading decay 
K+ 7r°e+z/. 

Consequently this decay can be reliably treated by the available field theoretical 
methods. 

The effective low-energy interaction mediating tt'^uj? is a fiavor-changing 

neutral current (FCNC). While being forbidden at the tree level, FCNC's can be 
induced through loops containing virtual heavy quarks. Therefore ti^vv 

• probes the standard model at the quantum level, thereby allowing for an 
indirect test of high energy scales through a low energy process, 

• is sensitive to the unknown top quark mass and 

• depends on the top quark couplings Vts and Vtd- 

Thus combined theoretical and experimental efforts may yield very desirable in- 
formation on parameters of the symmetry breaking sector in the standard model 
- the part of the theory least understood theoretically and containing the least 
known parameters such as the top mass and its CKM couplings. In particular it 
should be stressed that n^uu is probably the cleanest process to search 

for Vtd- 

The —>■ 71^ uu mode has not yet been observed experimentally. The published 



upper bound on the branching ratio is [11 



B{K+ 7r+z/z/)exp < 3.4 ■ 10" 
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A preliminary result from E787 at Brookhaven gives an improved value of 

B{K+ Tl^Uu)exp < 5 ■ 10"^ (2) 

which is still well above the theoretical expectation of O(10~^°), but a consider- 
able improvement is expected from upgrades planned for the coming years and 
designed to reach the 10~^° level |jl3[ . 



After these general introductory remarks let us now turn to a more detailed dis- 
cussion of the basic structure of the physics underlying ^ n^uu . The low 
energy effective hamiltonian relevant for the process K'^ n^uu may be written 
as follows 

G (y. 

'Heff = 2n ^ ^i'sVid X{Xi) {sd)v-A{i^i^)v-A (3) 

V 2 27r sm Uiy ^^^^^ 

where Xi = mf/M^r and X{x) is a monotonically increasing function of x to be 
given below. This hamiltonian consists of two parts: a top-quark and a charm- 
quark contribution, which originate to lowest order in the standard model from 
the Z-penguin- and box-diagrams shown in fig.l. By writing the hamiltonian in 
this form, the GIM mechanism is already incorporated and the up-quark contri- 
bution is understood to be subtracted from top- and charm-loops. In particular 
then X{x) has the property X(0) = 0. 

It is a characteristic feature of the Tt^vv decay amplitude that both top 

and charm contributions are of comparable size, since the smallness of X{xc) in 
comparison to X{xt) (due to rric <^ rrit) is compensated by the strong CKM sup- 
pression of the top contribution. By contrast this suppression is absent e.g. for 
the CP-violating mode Ki — > 7r°i/z/, where the branching ratio is determined by 
the imaginary part of the amplitude, or for rare B-decays such as -B ^ Xs^v and 
B l^l^ . Therefore the charm contribution is completely negligible in these 
latter processes. 

The short-distance electroweak loops of fig.l receive calculable radiative correc- 
tions through gluon exchange. Until recently these QCD corrections have been 
known only to leading logarithmic (LLA) accuracy |^, In this approxima- 
tion the charm contribution to the tt~^i>u amplitude is reduced by about 
35% through QCD effects. Since due to = 0(Mvi/) no large logs are present in 
the top contribution, this part remains uncorrected in LLA. As a consequence a 
non- negligible theoretical uncertainty remains due to the choice of the arbitrary 
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renormalization scale fi at which the top mass mt{fi) is defined. Although the ef- 
fect of varying /i around values of 0{mt) is formally of the neglected order 0{as), 
it amounts numerically to an uncertainty of 20%-30% for branching ratios. This 
issue has been discussed in detail in our previous papers ^ where we have 



calculated the complete 0{as) corrections to the top contribution for a class of 
rare decays governed by Z-penguin- and box-diagrams. This allowed us to give 
improved expressions, including the 0{as) corrections, for the branching ratios of 
those decays that are top-quark dominated {Kl tt^uu, B —>■ XsVv, B l^l') 
and to show that the scale ambiguities in the branching ratios are thus reduced 
from 0(25%) to typically 0(3%). 

The main objective of the present study is to perform a complete renormal- 
ization group analysis of the charm contribution to tt^uu in the next-to- 
leading logarithmic approximation (NLLA). Going beyond LLA is mandatory for 
consistency with the QCD corrected top-quark contribution. Furthermore it will 
decrease considerable theoretical uncertainties due to the choice of /i = 0(mc) 
present also in the charm sector. The remaining uncertainties will however be 
bigger than the ones in the top-sector due to the smaller scales involved at which 
as{fi) is larger and QCD perturbation theory less accurate. Combining the re- 
sult with the QCD corrected top contribution from |TB| will provide us with the 



final theoretical expression for B{K^ —>■ n^uu) which includes the short- distance 
QCD effects consistently at the next-to-leading-log level. 

At this stage we would like to recall the strong dependence oi B{K~^ n^uu) 
on the choice of as emphasized in particular by Dib [j^ and by Harris and 
Rosner ||^. These authors varying rric in the range 1.2GeV < rric < l.SGeV 
have found a 40-50% uncertainty in B{K^ n^uu) which is clearly disturbing. 
This large range for rric chosen in 0, H] is at first sight surprising because in a 
decay like K'^ n^uu the presence of the charm quark is only felt in the short 
distance loop and consequently it is evident that for the current quark mass 
value should be used. For the latter however one has mc{mc) = 1.3 ± O.lGel^ 



1^, On the other hand there is an ambiguity in the choice of the scale in 
rric and the use of mc{fi) with fi = 0{mc) in the calculation of B{K^ n^uu) 
is perfectly legitimate. With mc{m^ given above and IGeV < ^ < 3GeV one 
effectively ends up with a range for mc{fi) very similar to the one chosen for rric in 
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[0, |[. Consequently if one uses the leading order expressions as done in l]^, the 
uncertainty due to the treatment of rric is quite substantial. One of the important 
results of the present analysis is that the inclusion of next-to-leading corrections 
reduces considerably the uncertainty due to nic stressed in J^, |[. 

In order to further motivate the present calculation it is instructive to recall 
the relevant formulae for the top contribution |]15| and discuss their implications 
for the charm case. The function X{x) in is to 0{as) and without approxima- 
tion in m/My/ given by eqs. (^-([11]). In the case of charm (xc ~ 3 ■ 10~^ ^ 1) 
it is not necessary to keep all orders in m/Mw- We will restrict ourselves to 
keeping only terms linear in x, neglecting the order 0(a;^), which is an excellent 
approximation. Indeed in this case one can write the function X from @ with 
an error of less than 0.1% as (yU = m^; a = as/4:7r) 

N 3 , 1 / „ , 2 ^ , 23 + 27r2 

A X = — xinx x + a — 2xln x — 7xlnx x 

^ ^ 4 4 V 3 

It is instructive to take a closer look at the numerical values of the five terms in 
(^. For typical values = l.SGeV, Aqcd = 0.25GeV these five terms read 

X(x) = (16.32 - 0.66 - 13.05 + 5.54 - 1.37) ■ 10"^ = 6.78 ■ 10"^ (5) 

Several lessons can already be learned from this simple exercise: 

• One can see very clearly that straightforward perturbation theory is not 
reliable in this the 0{a) correction amounts to more than 50% of 
the lowest order result. 

• A renormalization group treatment is therefore mandatory. The LLA con- 
sists in summing the first and third terms in and similar terms 
0(xa" In*^"^^ x) to all orders in perturbation theory. 

• However the fourth term, neglected in LLA, is still quite sizeable. We 
can see explicitly that a next-to-leading log summation is very desirable: 
it resums all the terms of 0{xa"' In"' x), in particular the second and the 
fourth term above. 

• The non-leading mass term in the lowest order expression (— ;^x) is consis- 
tently taken into account within NLLA. 
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• The suppression of the charm contribution by leading-log QCD corrections 
is reflected in the negative sign of the 0{ax In^ x) correction term. The posi- 
tive next-to-leading log contribution indicates that this suppression tends to 
be weakened in NLLA. This will be confirmed by the exact result discussed 
in section 5. 

• The last term above is formally of the order 0{ax) and therefore neglected 
even in NLLA. Unfortunately it has a relatively large coefficient and gives 
roughly a 10% contribution. This number gives a hint on the remaining 
uncertainty to be expected. 

It is evident from this discussion that, given the high interest in the rare decay 
mode n^uu and considering the theoretical situation, it is worthwhile to 

perform the full next-to-leading order analysis in order to be able to reduce dis- 
turbing theoretical uncertainties. 

Our paper is organized as follows. Section 2 briefly collects some of the most 
important results. In sections 3 and 4 we will describe in some detail the formal- 
ism necessary for the two-loop renormalization group calculation in the charm 
sector for the Z-penguin- and the box-part respectively. In section 5 the results 
will be combined and the implications regarding the size of the corrections and 
the dependence on the renormalization scale will be investigated. In section 6 we 
combine both charm and top contributions and discuss phenomenological con- 
sequences for TT+i/P . In particular we address the determination of Vtd 
from this decay and we discuss the implications of our calculations for the uni- 
tarity triangle. The calculation of next-to-leading order corrections to the decay 
Kl — > A*"^/^"" is very similar to K'^ tt^vD and will be presented in section 7 for 
completeness. Unfortunately, due to insufficiently known long distance contribu- 
tions this decay is less useful for the determination of CKM parameters. We end 
our paper with a brief summary of the main results. 

2 The Main Formulae of this Paper 

Before entering the detailed discussion we find it useful to collect in this section 

the most important results of the present paper. 

First of all, the complete effective hamiltonian relevant for — > tt^uu can be 
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written as 



a 



eff 



72 27r sin^ Qw 



Y: {K:V,,X^jyL + V*Vt,Xixt))isd)v-A{m)v-A (6) 



=e,fj,,T 



The Ti^vv branching ratio obtained from is, for one single neutrino 

flavor, given by 



(7) 



Kf,27r2 sin^ Qw 

Here index / (/=e, /z, r) denotes the lepton flavor. The lepton mass dependence 
in the top-quark sector can be neglected. The function relevant for the top 

contribution, is to 0{as) and to all orders in m/Mw given by (a = Q;s/47r) 



with 



X^{x 



X{x) = Xo{x) + aXi{x) 
2 + X 3x — 6 



X 



+ 



1 — X (1 — x) 



Inx 



(9) 



and the QCD correction |T^, 

232; + — 4x^ x — llx^ + x^ + x'^ 



XAx) 



+ 



+ 8x 



3(1 -x)2 + 
+ Ax"^ + x^ — x'^ 

2(1 
(9Xo(x 



Inx 



m X — - —Loil — x) 



x] 



dx 



where x^ = /i^/M^ with /x = 0{mt) and 



^ , Int 
dt 



(10) 



:iii 



The /i-dependence in the last term in (p^Of ) cancels to the order considered the /x- 
dependence of the leading term Xq{x{^)). The corresponding expression for the 
charm sector is the function X]^^. It results from the RG calculation in NLLA 



and reads as follows: 



x[m) II 



xi 



32 



48 



NL 

24 



'.m 



696 



K- 



33 



+ I 1-ln^ (16ir+-8ir_ 



m 
/56248 
V 4375 



77 
1176244 
13125 



a(/x) 1875 

2302 3529184 

+ i^33 



6875 



48125 



81448 4563698 



6875 



144375 



(13) 
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where 



K = 


a{Mw) 




a{m) 




K_ = 


-12 

K— 


K33 



(1/2) _ x{m) M 



- In 




15212 



(14) 
(15) 



1875 



(1 - K 



625 



7500 



(16) 



Here K2 



rrir 



times omit the index / of X'j^j^. In (|T 



r = mf/m'^ and mi is the lepton mass. We will at 
- (|T6|) the scale /i = 0(mc). The two-loop 
expression for a{fi) is given in (P[). Again - to the considered order - the explicit 
ln(/i^/m^) terms in ( [13|) and ( [16|) cancel the /i-dependence of the leading terms. 
Expressing CKM elements through the Wolfenstein parameters A, A, q and t] 
with A = 0.22 and summing over the three neutrino flavors we finally find 



where 



Ti+vv) = 4.62- 10"i^v4*X^(xt 



(17) 



do 



1 + 



417 X'^L 



Formula (|T^ together with X{xi) and Xj^^ given in (|) and (|T^ respectively 
gives B{K^ n^uu) in the standard model with QCD effects taken into account 
consistently to next-to-leading order. The leading order expressions are obtained 
by replacing X{xt) Xo{xt) and X]^^ Xl with Xl found from (|13|) and 
(|TB|) by retaining there only the l/a(/i) terms. In LLA the one-loop expression 
should be used for a^. This amounts to setting /?i = in (0). A numerical 
analysis of ( |T7D and (p!8D together with related issues will be presented in section 



6. Numerical values of X'jsfi can be found in section 5. 



3 RG Calculation for the Z-Penguin Contribu- 
tion 

In the present section we will outline the basic formalism needed to perform the 
renormalization group analysis to next-to-leading order in the charm sector of 
the decay tt^vv. We start our discussion with the Z-penguin contribution 
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which is shghtly more comphcated than the box contribution. Having derived 
the former, which can be done in closed form, it will be staight forward to obtain 
the expression for the box contribution in the next section. 
The calculation consists of the following three steps: 

• In the first step the W and Z bosons being much heavier than the charm- 
quark are integrated out and the resulting effective hamiltonian is written 
down with matching conditions set up at /x = M\y- 

• The renormalization and mixing of the operators contributing to the sd — > 
Pu transition amplitude is calculated to two-loop order. The corresponding 
renormalization group equation is solved and the scale /i is evolved down 
to lower values. 

• At a scale n = 0{mc) the charm quark as an explicit degree of free- 
dom is removed from the effective theory. The only remaining operator 
is {sd)v-A{i^i^)v-A- Its coefficient yields the generalization in NLLA of the 
charm penguin contribution to X{xc). Since this operator has no anoma- 
lous dimension, no renormalization group running occurs below scales of 
0{mc) and the contributing QCD corrections are entirely perturbative. 

For the actual calculations we took external quarks massless and on-shell, which 
is the most convenient possibility. Since we are computing the short distance 
coefficients of an operator product expansion, the treatment of external lines 
is arbitrary and the final result does not depend on this choice. Furthermore 
we used dimensional regularization for both UV- and IR-divergences. We have 
employed Feynman gauge for the gluon and the W-boson field throughout the 
present work. 

3.1 Step 1 

After the W and Z bosons have been integrated out, the effective hamiltonian 
responsible for the Z-penguin contribution may be written as (Ac = V*gVcd) 
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where the operator basis is 



01 = -i j (fx T {{SiCj)v-A{Cjdi)v-A) {x) {{ckCk)v-A{vv)v-A) (0) - {c -> m} 

(20) 

02 = -i j d'^x T {{siCi)v-A{cjdj)v^A) {x) {{CkCk)v-A{^v)v-A) (0) - {c -> m} 

(21) 

0± = l(02±0i) (22) 

2 

Tfl 

Q = -^{sd)v-A{i^i^)v-A (23) 

For convenience we have introduced the bilocal structures Oi^2, which contain 
the relevant local 4-fermion operators in a form suitable for the evaluation of the 
sd vv transition amplitude. Furthermore we have switched to the diagonal 
basis 0±, anticipating the QCD renormalization and we have included a factor 
m? / g'^ (where m is the (charm-) quark mass and g the strong coupling constant) 
in the definition of Q in order to streamline the notation for the anomalous 
dimension matrices to be discussed below. 

As the first step we have to find the Wilson coefficients at to next-to-leading 
order. These can be written as follows {tF = (t>+, t>_, ^3)) 

v{Mw) = + a{Mw)v^^^ (24) 

^^°)^ = (1,1,0) (25) 
v^'^^ = iB+,B_,B,) (26) 



In the NDR-scheme {MS, anticommuting 75 in D 7^ 4 dimensions) one gets 

i?± = ±11^^ ^3 = 16 (27) 

where is the number of colors. The calculation of the scheme dependent next- 
to- leading order terms B± has been discussed in detail in |T^. The value of B3 
as well as the zero entry in ( P5| ) follows from matching the sd vv amplitude 
calculated in the effective theory (|19D onto the amplitude in the full theory to 0(a) 
at /i = Mvi/. To this end the renormalized matrix elements of the operators 0± 
are needed. They correspond to the diagram of fig. 2 and read in the NDR-scheme 

,2 \ 



= a(/i) ^7S(l- In ^)(Q) (28) 
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with 7j^3 given in (^Tj). Using (p8| ) and comparing the matrix element of (p!9D 
with the amphtude in the full theory 

Gf « . fx , 3 



Inx + -x \ {sd)v-A{i^i^)v-A (29) 



gives the value of quoted in (|27|). The amplitude A^^"^ follows simply from the 
evaluation of the Z°-penguin diagrams of fig.l in the limit x -C 1. 

3.2 Step 2 

In the next step we have to discuss the renormalization group evolution. The 
anomalous dimensions for 0±, to be denoted by 7-1- , are known to two-loop order 



and have been calculated in ||T8|, |I9| in various renormalization schemes. Here 
we use the NDR result of ||18[. The anomalous dimension of Q, denoted by 733, 
is also known to the same order. It depends only on the anomalous dimensions 
of the mass and the coupling constant as is evident from the definition (^31). 
What remains is the mixing of 0± to Q (7±3). The corresponding elements of 
the anomalous dimension matrix are obtained as usual from the divergent parts 
of the diagrams in figs. 2 and 3. The one-loop mixing has been considered in 
1^, ||. The two-loop case required for the present investigation is calculated 
in the present work and the results will be given below. In short, the anomalous 
dimension matrix has in the operator basis {0+, 0_, Q} the general form 

^ 7+ 7+3 ^ 
7= 7- 7-3 (30) 

V 733 / 

It governs the evolution of the Wilson coefficient functions by means of the renor- 
malization group equation 

V = 'j'^v (31) 



dln/i 

Note in particular that the evolution of v± is unaffected by the presence of Q. 
To proceed we briefly recall the general method for the solution of the renormal- 
ization group equation (^) in next-to-leading order, which we find useful also for 
the special case at hand. We follow [^]. 

The anomalous dimensions of the operators {0±, Q}, 7, and of the coupling con- 
stant, P{g)/g, are in second order of perturbation theory 

7(a) = 07^°) + 0^7(1) (32) 
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(33) 



We also need the anomalous dimension of the mass 

dm 



-1mm 



7m (a) = aimo + 



ml 



(34) 



(35) 



The solution of (|3T| ) can be expressed in terms of an evolution matrix ?7(/i, M^) 
as 

v{ii) = U{fi, Mw)v{Mw) (36) 
U{fi,Mw) = [l + a{fi)J]U'^'^\fi,Mw)[l-a{Mw)J] (37) 
U^^^ is the evolution matrix in leading logarithmic approximation: 

/ 



a(M, 



w 



a(/i) 



2/3o 



(3J 



D 



U^^\fi,Mw) = V 
where V diagonalizes 7^°)^ 

and 7*-'''' is the vector containing the diagonal elements of the diagonal matrix 

(0) 

Id ■ 

Defining now 



(39) 



G = V-'i^''^^V 
and a matrix S whose elements are 



(0) f^i 



(0) _ .,(0) 
j 



the matrix J is given by 

J = vsv-^ 

In our case the 3x3 matrix J has the structure 



(40) 



(41) 



(42) 



/ Jn \ 



J 



(43) 



J22 
\ J31 J32 J33 J 

similar to 7-^ in (0). 

We will now collect various expressions which enter the renormalization group 
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functions discussed above. 

For the QCD coupling in NLLA we use the form 

1 



a(/i) 



A2 



1 



A Inln^ 
PI In^ 



2 



(44) 



where A = A^^. The parameters of the /^-function and of the anomalous dimen- 



sion of the mass are (with N{f) the number of colors (flavors)) 

llA^-2/ ^ 34 _o 10 _ „ „ N"^-! 



3 ' 3 3 2A^ 

The nonvanishing elements of the anomalous dimension matrix read 



97 10 
7^1 = Cf ( 3Cf + - y / 



7^?= 2(7^^0 -/?o) 7^^=2(7^1 



(45) 
(46) 



7f = ±6^ 



(1) ^ /_2^ ± 57 ^ 19 ^ 4 ^ 



(47) 



7i1 = ±8(Ar±l) 



7^ = CFi±88N-4t 



The expressions 7*^^-* are in the NDR-scheme. The matrix V in (pOf) is 



/ 


1 





o\ 







1 





V 


1^31 


1^32 





31 



(0) 
7+3 



7+ 733 



7-3 



32 



(0) (0) 

7- -733 



Furthermore we define (see ([T^)) 



fCi = K 2.'30 



-,-,33 



(4J 



(49) 



(50) 



With all these formulae at hand we are now able to obtain the coefficients v at 
scales /i -C Mw by means of (0). We find 

v±ifi) = [1 + a(/i) J±] K± [1 - a(MH/)(J± - 5±)] (51) 

where = Ju, = J22 and 

t;3(/i) = 1^31 (i^+ - + V^2{K. - Kss) 

+a(/i) [(J3I + ^33^3l)i^+ + (^32 + J33V32)K^ - J33(1^31 + ^32)^^33] 

+a{Mw) [{B+ - Jii)V^iK+ + (fi_ - J22)V^2K- 

+ (53 - J31 - J32 + (^11 - 5+)V3i + (^22 - B_)V^2) K^z] (52) 
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3.3 Step 3 

The last step consists in integrating out the charm-quark. This is done by match- 
ing the effective theory in (|19D onto an effective theory without the operators 0± 
at a scale ^ = 0{mc). Below this scale the hamiltonian involves the single op- 
erator Q. The effects of the charm field being removed from the theory as an 
explicitly appearing degree of freedom are then still indirectly present in the mod- 
ified coefficient of Q. Effectively this step is accomplished by inserting (0±) of 
into dg). 



3.4 Final Result for the Z-Penguin 

Collecting the results of the three steps and using ([l9|) we obtain the final expres- 
sion in NLLA for the charm contribution to the Z-penguin part of the effective 
hamiltonian, suitable for scales below fi = 0{m^: 



'^eff,c 



a 



v/2 27rsin^eH' 



Ac Cnl {sd)v-A{j^Jy)v-A 



C 



NL 



x{p)_ 
32 



a{jj,) 



(53) 
(54) 



To the order considered the terms of 0(a) in f±(/i) have to be dropped. Catl is 
the generalization of the lowest order function Co(x) (see (|55|)) that includes the 
QCD corrections in NLLA. In fact 



X 3 

hni Catl = —Inx + -X = Co{x) 



(55) 



We will discuss this result in more detail later on. 

It is useful to express the running charm quark mass m(/i) in terms of m(m), 
which we use as an input parameter: 



X{fi) 



xim) 



a{fi) 




H 


a{m) 







(a(/i) — a{m)) 



(56) 



we then obtain formula ( |13|) quoted in section 



Setting iV = 3 and / = 4 in (I 
2. 

We have performed the renormalization group evolution in an effective four-flavor 
theory from scales of 0{My/) down to 0{m^ and neglected the effects of a b- 
quark threshold and a corresponding five-flavor theory above m?,. This is an 
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excellent approximation which allows one to write down a fairly compact final 
result in closed form. We have checked that the inclusion of a b-quark threshold 
would change the values for the charm-quark function X^l by not more than 



3.5 Renormalization Scheme Independence 

To conclude this section we turn to the question of renormalization scheme de- 
pendences and their cancellation. This issue is well known and has already been 
discussed elsewhere (see e.g. [^). However since the application of the next- 
to-leading renormalization group formalism to the case considered here is new 
and the explicit expression oi Cnl clearly involves scheme dependences whose 
cancellation may not be entirely obvious from (0), we find it useful to give a 
short account of this subject in explicit terms. 

In essence the scheme independence of a matrix element of the effective hamilto- 
nian (|19D is trivial since it is proportional to 



which is manifestly scheme independent. Different choices of the scheme just shift 
finite 0(a)-terms between the coefficients and the operator matrix elements- 
this freedom is in fact the reason why the unphysical coefficients are scheme 
dependent. More explicitly these quantities are in different schemes related by 



0.1% 



(d^(/i))t;(/i) = (i;+0+ + v^O^ + v^Q)\ 



(57) 



{d)' = {l-as){d) 



(1 + as^)v 



(58) 



where s is a constant 3x3 matrix of the form 



/ 




s = 



s 



(59) 



V ss J 



Since 



(d^(/i))t;(/i) = {d^ifx))U{fx,Mw)v{Mw) 



(60) 



implies 



U'ifi, Mw) = (1 + a(/x)s^)[/(/i, Mw)il - aiMw)s^) 



(61) 



we have from ( p7|) 



J' = J + s'^ 



(62) 
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If we use the relation between J and 7'-^'' (from (^D|) - (^)) 



J + 



+ £507^"^' (63) 



2/5o 2(3', 



we obtain the transformation property of the anomalous dimensions under a 
change in the renormalization scheme: 

^(1)' = ^(1) _ [s, 7^°)] - 2l3os 7^°^' = 7^°^ (64) 

Given these general formulae one can easily clarify the cancellation of scheme 
dependences in all particular cases or convert scheme dependent quantities from 
one scheme to another, if desired. To be more specific, let us consider expression 
(p2D. The first line being a function of one-loop anomalous dimensions and the 
coupling constant only is obviously scheme independent. The coefficient of a{Mw) 
is scheme independent since the scheme dependences in J and Bi cancel in the 



combinations 5+ — Jn, — J22 and i?3 — J31 — J32 ((^), ( ^8]) and (|6^)). The 



scheme dependences in the coefficient of a(/i) remain and are responsible for the 



scheme dependence of v^{jji). We have, using (|52D , (|62D and {Q)' = (1 — as^){Q) 
v'siQY = v,{Q) + a(/i)[s+3ir+ + s_;K_]{Q) (65) 
On the other hand, up to terms of the neglected order (note that {0±) = 0(a)) 
v'+iO+y + v'_{0.y = v+{0+) + v.{0.) - a{fi)[s+3K+ + s.,K_]{Q) (66) 



Eq. (pBD reflects the arbitrariness one has in the subtraction of finite constants 
when defining the renormalized matrix elements {0±). The scheme dependences 
in ( |65D and (|66|) cancel in the sum, which essentially represents the function C^l. 
The cancellation of the /^-dependence will be discussed in section 5. 



4 RG Calculation for the Box Contribution 

After the detailed exposition of the Z-penguin contribution in the preceding sec- 
tion it is straightforward to repeat the analysis in the analogous, but somewhat 
simpler case of the box. We will briefiy summarize the essential steps. 
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4.1 Step 1 

When the W boson is integrated out, the hamiltonian relevant for the box diagram 
is given by 

V2 27rsin^9i4/ \ J 
= -ijd^xT ((sc)v-A(^>/)y-A) {x) [ilu)v-Aicd)v-A) (0) - {c ^ m} (68) 

2 

Q = -^{sd)v-A{i^i^)v-A (69) 

In order to find the initial conditions for the coefficients Cj we need the renormal- 
ized matrix element of the operator O. Calculating the diagram of fig.4 we find 
in the NDR-scheme 

(0) = a(/.)16(ln^ + ^ + [^)(g) (70) 

Inserting this result into ( |67D and comparing with the amplitude in the full theory 
given by 

= -T2^^^'' (-1- + - (l + ^)) (■5<')v-.«^W.-. (71) 
implies the initial values for the Wilson coefficient function in the NDR-scheme 
c'iMw) = (c.iMw), C2iMw)) = (1, 0) + aiMw)iO, B^) = -36 (72) 



4.2 Step 2 

The anomalous dimension matrix in the basis {O, Q} has the form 

V 722 / 

with the perturbative expansion 

7 = 07^°) + 0^7(1) (74) 

The anomalous dimensions '^^21 Ivii describing the RG mixing of operator O 
with Q follow from the divergent parts of the diagrams in fig.4 and fig. 5 respec- 
tively. 7i2^ has been calculated in ||, |^. 7(2'' is new. The solution of the 
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renormalization group equation in NLLA is obtained via the general method de- 
scribed in section 3. 

The matrix V needed to diagonahze 7*^°)-^ 

V-^,mv = 1' ° I 1 (75) 



is given by 



7g' 



V2I 1 / 722 



Furthermore, the matrix J defined in (^2]) can in this case be written as 



The nonvanishing elements of the anomalous dimension matrix are 

72? = 2(7mo - Po) I22 = 2 (7ml - A) 



(7J 



7i? = -32 7i2^ = 80C^ 



Defining K2 = ii"''22'' /2/3o then find for the solution of the renormalization group 
equation 

-^c = 7^c- (79) 
ctln/i 

ci(/i) = 1 (80) 
c^{^) = 1^21(1 - K2) + aifi) [J21 + ^21^22(1 - K2)] + a{Mw){B2 - J2i)i^2 (81) 



4.3 Step 3 

The last step consists of integrating out the charm quark. This adds an additional 



contribution to the coefficient of the operator Q which in view of (|80|) is simply 
given by the renormalized matrix element of the operator O at a scale ^ = O (m) , 
already quoted in ([TOl). 

4.4 Final Result for the Box Contribution 



Collecting all this and using (pTf ) we finally obtain the effective hamiltonian in- 
duced from box diagrams with charm, including the NLL QCD effects between 
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the scales Mw and /i = 0(mc): 

^S,c = ."n (.-rf)v-^(z7,z.0y-A (82) 

V2 27rsm HvK 



p(i/2) _ a;(/i) 



16 fln4 + ? + ^U'^^^^ 



3) 



4 1 — r J a{n) 

(^2]) is written here for one z/-flavor. The index (1/2) refers to the weak isospin 
of the final state leptons. 

If we set = 3, / = 4 and use ( p6D we arrive at the expression given in (plGl). 
We would like to add a few important remarks. 

• As usual we are working in the NDR-scheme; all scheme dependent quanti- 
ties that we have given above (one-loop initial values and two-loop anoma- 
lous dimensions) refer to this scheme. However in the case of the box con- 
tribution there is an additional scheme dependence involved which arises 
from the treatment of the particular Dirac structure. In our calculation we 
have used the following projection 

7V7'^(1 - 75) ® 7^.7.7.(1 - 75) ^ 4(4 - e)ni - 75) ® 7,^(1 - 75) (84) 

A different treatment of this expression would lead to different values for 
B2 and 712'' and would correspond to a different renormalization scheme. 
This "projection scheme dependence" of course cancels in the final result 
provided coefficient and anomalous dimension are evaluated consistently. 

• The proof of the cancellation of all scheme dependences (e.g. the projection 
scheme dependence mentioned above) proceeds in the same way as discussed 
at the end of section 3. 

• The operator O contains quark fields only in the form of local currents 
and therefore has no anomalous dimension. This feature is reflected in the 
structure of the anomalous dimension matrix in (|73D . As a consequence the 
coefficient Ci of O is /^-independent (^0]). 



In contrast to the case of the Z-penguin the box contribution involves an 
additional mass parameter mi, the mass of the internal charged lepton in 
the box diagram. While mi can be neglected compared to the charm mass 
(r ^ 0) in the case of the electron or the muon, it is certainly non-negligible 
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in the case of the tau. The effect of a nonvanishing lepton mass of the order 
of the charm quark mass is represented by the r-dependent term in the 
expressions for -B^{^'* in (|83D , (|16]). A comparison with (^) shows that the 
dependence on m; arises entirely from the matrix element of the operator 
O, evaluated at the scale fi = 0{mc). Note in particular that the initial 



values for the Wilson coefficients (^) and the anomalous dimensions (|78D 
do not depend on mi. This is intuitively clear since the Wilson coefficient 
depends only on the short distance structure of the effective interaction and 
must therefore be independent of the low energy scale mi. 
In this context the GIM mechanism plays an essential role. The contracted 
box diagrams with charm (figs. 4, 5) obviously have a quadratic divergence. 
However, what really matters for the matrix element of O is the difference 
between the charm- and the up-quark contribution (see (|68D). As a conse- 
quence the matrix element (O) is only logarithmically divergent. It follows 
then that the difference between the cases mi = and m; 7^ is finite and 
hence does not affect the anomalous dimensions, as it must be in a mass 
independent scheme like minimal subtraction. To conclude, the effect of 
the GIM mechanism is nontrivial in the case of the box with nonvanishing 
lepton mass, whereas it is not explicitly visible for mi = (and in the case 
of the Z-penguin) since in these latter cases the up-quark contribution is 
anyway vanishing in dimensional regularization. 

5 Discussion of the Charm Contribution 
to Next-To-Leading Order 

In sections 3 and 4 we have calculated the Z-penguin and box-parts of the charm 
contribution to the effective hamiltonian for —>■ n^uu decay in NLLA. The 
explicit results for the functions Cnl and -B^^,^^ are given in ([T3|) and ([T6|) respec- 
tively. Combining these functions, which depend on the W-boson gauge, into the 
gauge independent combination 

Xnl = Cnl — ^B^NL^ (85) 
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we can next write down the final result for the charm contribution to the effective 
hamiltonian 

G F 01 

T-ieff,c = —7=- — ■ 2n isd)v-A{i^i^)v-A (86) 

We proceed to discuss various aspects of this result in what follows. 



Let us first compare the charm- with the top-contribution as obtained in |T^ . In 
the case of the top the internal quark mass is large and comparable to M\y so 
that no large logarithms appear. The charm mass on the other hand is very much 
smaller than Mw and logarithmic terms become important. For these reasons the 
two cases, regarding the theoretical methods required, are quite complementary 
to each other: For the top-part a consideration of QCD effects in straightforward 
perturbation theory to 0{as) is appropriate, but all orders have to be kept in the 
mass ratio rrit/Mw- On the contrary the charm contribution allows a restriction 
to the order 0{xc) in the mass ratio but requires summation of leading and next- 
to-leading logarithms to all orders in QCD. From this consideration it follows 
that one can check the top contribution against the charm contribution by taking 
the limit x <^ 1 in the former case and simultaneously expanding the RG result 
to 0{a) in the latter. 

Recall that the function X in (H) consists of a Z-penguin- and a box-part C and 
B, respectively 

X{x) = C{x) - 4B{x, +1/2) (87) 



Keeping only terms linear in x in the exact results of |TH] we have {mi = 0) 



^/ N . 1 1 3 / , 2 23 , 29 + 27r2 \ 

C{x) = -xmx-\ — x + a 2xm x -\ xlnx-\ x] 88 

^ ^ 4 4 V 3 3 J ^ ' 

/ 44 52 + 47r2 \ 
4i?(x, +1/2) = xlnx + x + a ( 4xln x + —x\nx-\ ^ xj (89) 

Expanding now the RG-functions Cnl and fij^^f ^ (for r = 0) in (0) and (|16D to 
0(a), thereby "switching off' the renormalization group summations, we obtain 
in fact the first four terms in ( |88| ) and ( pQ]) respectively. Let us mention a few 
interesting points. 



The last terms in (|88|) , (|89|) are of 0{ax) and go beyond the NLLA; therefore 
they can of course not be obtained from Cnl and B^l- 
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A comparison of Bnl, expanded to 0(a), with (^) may also be used to 
determine the two- loop anomalous dimension '~f[^2 ■ Once 7^2"* is known, the 
full next-to-leading order expression for the box can be obtained. This 
observation allows to circumvent the explicit evaluation of the two-loop 
diagrams in the effective theory (fig. 5) - the information on 7I2'' is already 



contained in the calculation of B{x, +1/2) in [jTS 



• By contrast, the same exercise for Cnl shows that (^) only yields the 
value for (7,^3 +7i"'^3). To get both elements 7^3 separately it is necessary to 
investigate the calculation in the effective theory (fig. 3). Still this feature 
can be used as a cross-check. 

An important issue is the independence of the physical amplitude on the renormal- 
ization scale fi. Since in our case the operator {sd)v~A{i^^)v-A has no anomalous 
dimension, also the coefficient functions C, B and X are //-independent. How- 
ever, because we are working in perturbation theory, the //-independence is valid 
only up to terms of the neglected order. As a side remark we mention that the 
renormalization scale /t appearing in our final expressions Cnl snd Bf^i^ can be 
viewed also as the scale at which the charm quark is removed from the effective 
theory. This interpretation is clear from step 3 of the derivation given in section 
3. Although this scale should be of the same order as the charm quark mass, 
it is not necessary that /i equals rrtc exactly. Therefore we have displayed the 
/i-dependence explicitly in our formulae and have not set /i = nic. This allows us 
in particular to investigate the residual /i-dependence. To do so consider first the 
leading-log parts of Cnl and B^l (in O) and (|^)) and their variation with /i: 

d ( x(/i) 



[F3l(ir+-ir33) + V^32(^--^33)] 
X 

32 



c?ln/i y32a(/i) 

[7?]^+ + 7i1i^- +0{ax) (90) 



d ( x(/i) \ X (0) 

— — — V2I [Ko - 1 = -— 7l9 



We can see that the residual /t-dependence for the leading-log result is of the 
order 0{x). Formally terms of this order are neglected in LLA and the results 
are /i-independent to the order considered in this approximation. Numerically 
however the residual /i-dependence is quite pronounced as we will see below. 
Furthermore we note again that treating QCD corrections in LLA and keeping 
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the non- logarithmic 0{x) terms in Xq{x) would strictly speaking be inconsistent; 
a variation of the scale fi, which is not fixed in LLA, would give rise to terms 
of the same order. Similarly the non-negligible contribution from the tau-lepton 
mass term is a next-to-leading order effect (see (p3D). 

All these considerations underline once more that the leading-log approximation 
is unsatisfactory. An improvement can only be achieved by going beyond the 
leading order. In particular in NLLA the /^-dependence of 0{x) in the leading 
order terms of Cnl and Bj^l given in ( pOD and ( pT]) is cancelled by the explicitly 
//-dependent terms appearing at next-to-leading order, as is evident from (Q) 
and (|83D . As a result the remaining /{-dependence is 0{ax) instead of 0{x). 
What this amounts to numerically is illustrated in fig. 6 where we have plotted 
the charm-function Xjyi (for mi = 0) compared to the leading-log result and 
the case without QCD for values of /i between 1 and 3 GeV. We observe the 
following features: 

• The residual slope of X^vl is indeed considerably smaller in comparison to 
Xl, which exhibits a quite substantial dependence on the unphysical scale 
/i. The variation of X (defined as {X{lGeV) -X{3GeV))/X{m,)) is 18.8% 
in NLLA compared to 52.7% in LLA. 

• As already anticipated in the Introduction, the suppression of the uncor- 
rected function through QCD effects is somewhat less pronounced in NLLA. 

• The next-to-leading effects amount to a ~ 10% correction relative to Xl 
at /i = rric- However the size of this correction strongly depends on fi due 
to the scale ambiguity of the result. This means that the question of how 
large the next-to-leading effects compared to the LLA really are cannot be 
answered uniquely. Therefore we focus our attention on the more important 
issue of the /i-dependence and its reduction in NLLA. 

Whereas it is meaningful to use A^g^ in Xj^l one has no control on the exact 
definition of the QCD scale in leading order (Alo)- Alo must be of the same 
order of magnitude as Ajjg but is undetermined otherwise, which constitutes 
a further problem of the LLA. For definiteness of the comparison we have set 
Alo = ^JTs ~ 0.25GeV in fig.6. Finally in table 1 we give the numerical values 
for Xffj^ and X'^j^ for fi = rric and several values of A;^ and mdrric). The 
reduction of Xnl through the effect of rrir is clearly visible. 
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^MS' [GeV] 


1.2 


1.3 


1.4 


1.2 


1.3 


1.4 


0.15 


10.24 


11.95 


13.78 


6.99 


8.40 


9.92 


0.25 


9.51 


11.18 


12.97 


6.26 


7.63 


9.11 


0.35 


8.74 


10.39 


12.15 


5.49 


6.83 


8.29 



Table 1: The function Xnl for various and rric- 

6 Phenomenological Implications for 

TT^uu, Vtd and the Unitarity Triangle 

6.1 General Expressions 

In the present section we will finally discuss the expression for the branching ratio 
of TT'^uu , including the complete 0{as) corrections to the top contribution 

and the renormalization group result in NLLA for the charm part. As applications 
we will discuss consequences for a possible determination of | VJ^ I from a measured 
B{K~^ —>■ Ti^uu) and the impact of our results on the unitarity triangle. In our 
presentation we follow chapter 6.2. of [|I|. There a detailed investigation of the 
Ti^uu phenomenology has been performed with QCD corrections treated 
in LLA. While the discussion in [Q] concentrated on the various uncertainties 
arising from standard model parameters like CKM elements and quark masses, 
we will put emphasis on the theoretical uncertainties due to residual dependences 
on the renormalization scale. 

Taking (|^), summing over the three neutrino fiavors and expressing CKM vari- 
ables in terms of the Wolfenstein parameters A, A, g and rj, we obtain for the 
branching ratio 



27r2 sin^ 9 



w 



where 



^0 = 1 + 



[ 2 ^NL 



(92) 



(93) 
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Using the numerical values for various parameters as collected in ( |1U1| ) and ( |102D 



we find the expressions quoted in section 2. Here and in the following we have 
used the QCD corrected top contribution X{xt) as it stands and made no effort 
to expand functions of X{xt) to first order in a^. For renormalization scales not 
too far from rrit this is anyway irrelevant numerically. 
Eg . (1921) defines a circle in the {g,ri) plane with center at {go,0), where 

\2 ve 1 ^ VT 



and with the radius squared 



1 

r„ — 



2 \2 



X^A^X^ixt] 
Here we have defined 



sin"" ewB{K+ ^TT+uu) 2 / A 
3a^B{K+ ^ TT^e+z/) 9 [^^ 2 ^^^^^ 



(95) 



XT 111 T 

r — 1 



The restriction to leading logarithmic accuracy in the above formulae is made by 
replacing X{xt) -^0(2;*) and X^j^j^ X^. The difference SX'^j^ should be set 
to zero in LLA for consistency as discussed in section 5. 

We recall a few important definitions related to the unitarity triangle in the {g, if) 
plane. First we introduce |jl|] 

i?, = = JY!^ (97) 

To a very good approximation these can be expressed in terms of g and 1] as 



Rb = ^/g^ + V^ Rt = ^{l-g? + V^ (98) 

The circle defined in the (^, r]) plane by the value for B{K^ n^uu) ((|S^) and 
(p5|) ) intersects for the allowed range of parameters with the circle determined by 
a value for Ri,. This allows to determine Rt. One finds 

R^^ = l + Rl + d^-g^ (99) 

go 

6.2 The Analysis of 1^^^! 
The modulus of Vtd is obtained as 

\Vtd\ = AX'Rt (100) 
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mt/GeV 


mJGeV 




A 


Rb 


B{K+ ^ 7r+z/p)/10-i° 




150 


1.3 


0.25 


0.89 


0.45 


1.0 


I 


±50 


±0.1 


±0.10 


±0.08 


±0.14 


±0.5 


II 


±5 


±0.05 


±0.05 


±0.04 


±0.05 


±0.1 



Table 2: Collection of input parameters with two different choices for the un- 
certainties. 



We note that besides A (or |Kb|) and rrit the value of -Rf, (or iKtfe/K^I) is needed 
together with the branching ratio in order to extract \Vtd\ from ti^vv 
decay. 

Before discussing numerical examples we will now compile the input parameters 
necessary for the analysis. We take 

Mw = SOGeV rrir = 1.78^6^ A = 0.22 (101) 

a = 1/128 sin^ = 0.23 5(ir+ ^ 7r°e+z/) = 4.8 ■ 10"^ (102) 

The remaining parameters are given in table 2 together with two different sets 
of uncertainties. Range I corresponds to the present status with the exception of 
B{K^ TT'^uu), where we have just chosen a typical number. For later discussion 
we have also displayed range II illustrating some possible future scenario. The 
values of A and Rf, correspond to = 0.043 ±0.004 and iKb/Kfel = 0.10 ±0.03 
(range I) , which are in the ball park of various determinations. The mass param- 
eters rrit and rric in table 2 stand for mt{mt) and mc{mc). 

The various uncertainties in all these parameters result in a rather large range of 
predicted values for B{K^ tt^uu), typically (1 — 6) ■ 10^^°. Similar uncertain- 
ties exist for the determination of \Vtd\ from the branching ratio |jl], 0. 
However, these uncertainties from input parameters should be distinguished from 
the uncertainties arising due to the /i-dependences of the physical amplitudes, 
which have their origin in the truncation of the perturbation series. The param- 
eters could in principle be considered as determined from somewhere else and in 
the ideal case known without errors. However even in this case there would still 
remain the /i-scale ambiguity which constitutes the intrinsic uncertainty of the 
theoretical prediction. An improvement in this issue would require going even 
beyond NLLA. 

For these reasons we will focus our attention on the scale ambiguity and compare 
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LLA and NLLA in this respect. As an illustration of this point let us consider 
the determination of \Vtd\ via (|100|) . We fix all input parameters at their cen- 
tral values but allow a variation of the renormalization scales within the ranges 
IGeV < fic < ^GeV and lOOGeV < ^it < SOOGeV" in the charm- and the top- 
sector, respectively. The variation of the unphysical scales /ic and ^it is of course 
not unique. To avoid large logarithms and /i^ should however not deviate too 
much from tjIc and respectively. Our ranges of fJ't, chosen for definiteness, 
satisfy this requirement. The ensuing variations in predictions of physical quan- 
tities are to be understood as indications for typical theoretical uncertainties to 
be expected due to the truncation of the perturbation series. 
We find 

8.58- 10"^ < l^tdl < 11.35- 10-3 LLA (103) 
9.51 ■ 10-3 < \Vtd\ < 10.22 ■ 10-3 NLLA (104) 
The central values at fic = ttLc and fit = are 

= 9.49- 10-3 LLA (105) 

|V^d| = 9.68- 10-3 NLLA (106) 

One observes that including the full next-to-leading order corrections reduces the 
total variation from 29% (LLA) to 7% (NLLA) in the present example. The main 
bulk of this ambiguity stems from the charm sector. Keeping /ic = irtc fixed and 
varying only /if, the ambiguities would shrink to 10% (LLA) and 1% (NLLA). 
These purely theoretical uncertainties should be compared to those coming from 
the input parameters. To get an idea of their magnitude we show in table 3 the 
relative change of \Vtd\ when one parameter is varied in the two ranges given in 
table 2 while all other parameters are kept fixed. The entries in the first row 
indicate the varying parameter. We note that the dependence on I^jjg and Rh is 
quite small. It is clear that at the present stage with B{K^ n^uu) being un- 
known and with a large range of possible top quark masses, the theoretical errors 
related to scale ambiguities seem rather marginal. However the determination of 
relevant parameters should improve in the future. Once the corresponding preci- 
sion will have attained the level exemplified by our choice of parameter range II in 
table 1, the resulting uncertainties will start to be comparable to the theoretical 
uncertainty of about 7%. Clearly, in such a situation the gain in accuracy of the 
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^\Vtd\/\Vtd\ 


mt 


rric 




A 


Rb 


B{K+ -K+uu) 


I 


78.5% 


11.6% 


5.7% 


16.8% 


3.4% 


70.8% 


II 


7.2% 


5.8% 


2.8% 


8.4% 


1.2% 


13.3% 



Table 3: Relative sensitivity of \Vtd\ to variations of input parameters according 
to table 2. 

theoretical prediction by a factor of more than 4 compared to the LLA will be 
very important. In order to show this more explicitly we plot in fig. 8 \Vtd\ as a 
function of mt{mt), keeping all other input parameters at their central values but 
varying /ic and fit according to IGeV < fic < SGeV and |mt < fit ^ 2mt. The 
broad and the narrow band correspond to LLA and NLLA respectively. This fig- 
ure illustrates very clearly that the next-to-leading QCD analysis is necessary to 
really exploit the clean nature of tt'^uu in future constraints on important 

standard model parameters. 

To summarize: Whereas in LLA the scale ambiguity would be the main uncer- 
tainty in \Vtd\ provided mt and B{K^ n^uu) have been accurately measured, 
the main uncertainty after our calculation of next-to-leading order corrections 
comes from the value of Vet or A. 

6.3 The Impact on the Unitarity Triangle 

It should be emphasized that the measurement of B{K^ n^uu) will not only 
determine \Vtd\ but also constrain the shape of the unitarity triangle. Indeed the 
parameters g and r] can be determined from B{K^ rr^uu) once Rb, \Vcb\ and 
rrit are known. We find 

V = -JRI - t!^ (108) 

with Qo and Tq given by ( pif ) and (^51) . r] is bound to be positive when the 
experimental result for the CP violating parameter e is taken into account [jl|. 
Here we would like to discuss the sensitivity of Tq, Qo, Q and rj to the choice of 
renormalization scales. We first note following [jl| that the second term in (^) is 
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negligible. Consequently 



° A2X(xOV 4.62-10-11 ^ ^ 

is entirely determined by the top contribution to — n^uu. Now as we stressed 
in [jl4|, and in discussing \Vtd\ above, the scale dependence in X{xt) is after 
the inclusion of next-to-leading order corrections reduced from 0(10%) to 0(1%). 
Consequently the theoretical calculation of r^, is fully under control once A, rrit 
and B{K^ tt^uu) have been precisely measured. 

In the absence of charm contributions one would have = 1- Thus the de- 
parture of Qo from unity measures the impact of the charm contribution on the 
determination of the unitarity triangle. Since in LLA there was a considerable 
scale ambiguity in the charm contribution we expect this to show up also in the 
determination of Qo, g and rj. To summarize this discussion we have displayed in 
table 4 the /i-dependences of Tq, Qo, g, r], \Vtd\ and B{K^ vr+z/i/). Again we 
take the central values for all input parameters but vary /ic and fit in the ranges 
IGeV" < fic < 3GeV, lOOGeV" < fit < SOOGeV". For each quantity we show 
the resulting minimal and maximal value, in NLLA as well as in LLA, together 
with the percentage of this variation refering to the standard value (/ic = rric 
and fit = rrit). For the evaluation of B{K^ n^uu) one has to fix g and 
7]. We have chosen g = 0.08, rj = 0.44 which correspond to a central value of 
B{K+ TT+z/z/) = 10-1° in NLLA. 

As anticipated the scale uncertainties present in LLA are considerably reduced in 
NLLA. This is in particular seen in the case of g where the uncertainty Ag = 0.31 
present in LLA is reduced to Ag = 0.08 after the inclusion of next-to-leading or- 
der corrections. The very small ambiguity in rj results partly from the small value 
of g chosen in this example, but generally rj is less sensitive to fi than g. 
In fig. 9 we show the position of the point {g,ri) which determines the unitarity 
triangle. To this end we have fixed all parameters at their central values except 
for Rh for which we have chosen three representative numbers, i?b=0.31, 0.45, 
0.59. The full and the reduced ranges represent LLA and NLLA respectively. 
The impact of our calculations on the accuracy of determining the unitarity tri- 
angle is quite impressive. We note that the reduction in the uncertainty in g may 
allow to resolve the two-fold ambiguity {g > 0, g < 0) present in the analysis of 
the CP-violating e-parameter. 
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NLLA (LLA) 


ro 


Qo 


Q 


V 




1.39 (1.37) 


1.39 (1.39) 


0.08 (0.10) 


0.44 (0.44) 


mill 


1.38 (1.31) 


1.34 (1.24) 


0.02 (-0.12) 


0.44 (0.41) 


max 


1.40 (1.45) 


1.41 (1.49) 


0.10 (0.19) 


0.45 (0.45) 




1.0% (9.8%) 


5.3% (17.6%) 




2.3% (9.7%) 



NLLA (LLA) 


\Vm\/io-' 


B{K+ 7r+i/z>)/10-i° 




9.68 (9.49) 


1.00 (1.03) 


mill 


9.51 (8.58) 


0.92 (0.76) 


max 


10.22 (11.35) 


1.03 (1.19) 




7.4% (29.2%) 


10.9% (41.0%) 



Table 4: Various quantities relevant for the phenomenology of vr+z/z/ and 

their theoretical uncertainties in LLA and in NLLA. 

We note also that the 40% uncertainty in B{K^ vr+z/z/) present in LLA is 
reduced to 11% in NLLA. This is also seen in fig. 10 where B{K'^ —* vr+z/z/) is 
given as a function of rrit for central values of the parameters. The meaning of 
the curves is as in fig. 8. 

7 Short Distance Contribution to Kl 
Beyond Leading Logs 

7.1 Renormalization Group Analysis 

The branching ratio B{Kl — > f^^f^~) has already been measured in experiment 
pl[l . However, since this decay receives important long distance contributions, a 



complete theoretical treatment is more difficult than in the case of tt^vv. 



Recent discussions of this issue have been given in ||22|, ^ . Yet the pure short 
distance part of Kl f^^fJ'~ is on the same footing as —>■ tt^uu and can be 
calculated by the same methods. Since the computation is very similar in both 
cases, we will restrict ourselves to mainly quote the final results for -—>■ fj^^fi~, 
indicating the most important differences to the case of tt^uu. 
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The effective Hamiltonian inducing the short distance contribution to Kl 
can be written as 

G p 



n 



a 



eff 



2n sin^ 6 



w 



(KtKdlVL + VZVt,Y{xt)) {sd)v-A{m)v-A (110) 



The resulting branching ratio is 



The function Y{x) is given by 



Y{x) = Yo{x)+aYi{x) 



where 



Yn(x 



A - X 



+ 



3x 



\nx 



and 



Yi(x) 



1 — X (1 — x)^ 
Ax + 16a:^ + 4x — lOx^ — x^ — 



111) 



(112) 



;ii3) 



ad -x) 



(1 -X)2 



Inx 



+ 



2x — 14x^ + x^ — x^ , 2 . 2x + x'' 



+ 8x 



2(1 -x)3 
dYoix) 

OX 



In X + 



1-x)- 



Xofl-X) 



The RG expression Y^v^ representing the charm contribution reads 



Y, 



NL — '-^NL 



Cnt, — B 



(-1/2) 
NL 



(114) 



;ii5) 



where Cat^ is the Z-penguin part given in (|TBp and B^l^"^^ is the box contribution 
in the charm sector, relevant for the case of final state leptons with weak isospin 
T3 = -1/2. We find 



B 



(-1/2) 
NL 



x(m)^|| 



3(1 - Ko 



a{i2) 1875 



jj,^ 329 15212 _ 30581 



- In ^ h 

^2 12 625 



-K2 + 



7500 



KKo 



(116) 



Note the simple relation to B^^^^^ in (|T6[) (for r = 0) 



B 



(-1/2) _ ^(1/2) x{m^ M 



NL 



'NL 



Kc'' {KK2 - 1) 



(117) 



30 



For the Z-penguin contribution the change from neutrinos to muons as exter- 
nal leptons is trivial and the corresponding function Cnl is the same in both 
cases. The box-part on the other hand is slightly different for Kl —>■ fi~^fi~ and 
— > vr+z/i/ since the lepton line in the box diagram is reversed. As a first con- 
sequence no internal lepton mass appears in the case of Kl —>■ ^~ . The renor- 
malization group calculation is however very similar to the one for —>■ tt^vv. 
We comment briefly on the most important differences. Starting point are again 
eqs. (p7|)-(|69D but with the operator Q now defined as 

,2 



Q 



m 



9' 



V-A 



118) 



Using the NDR-scheme and the projection 

r^^Yi^ - 75) ® lulalf.{l - 75) ^ 4(1 - 2^)7^(1 - 75) ® 7m(1 - 75) (119) 
we have instead of ( [7DD 



(0) = -a(/i) 41n^-2 (g) 



(120) 



For 5^2^^'* o'^s then finds 



B 



(-1/2) 
NL 



x(/i) 



-41n4 + 2 + ^^ 



1211 



The Wilson coefficient C2(/i) is as given in (|8lD , however with different values for 
B2 and 712: 



B2 = 2 7i? = 8 7I2' = 160C^ 



(1) 



122) 



Putting all this together, using ( ^61) and setting = 3, / = 4 one finally obtains 
(i™. 



7.2 Discussion of the charm contribution 

We now turn to the discussion of theoretical uncertainties. In the top-quark 
sector the residual scale dependence is almost completely eliminated when the 
0{as) corrections are included ( |112| ), just as in the case of n^uu. The 

charm contribution is again more problematic in this respect. In addition it 
exhibits a somewhat different structure than the charm sector in tt^uu. 



31 



To see this, recall that the function Y |r3| reads in terms of the Z-penguin- and 
box-contribution 

= C(x) -5(x,-l/2) (123) 
where, to linear order in x, C is given in ( ^8|) and 

—1/2) = -xlnx + -X + a ^xln^x + -^xlnx H — — (124) 

It follows that the terms of O(xlnx) and 0{ax\n.x) both cancel in 

1 / 4 + 7r2 \ 
Y{x) = -^x + a\x\v? X -\ ^ — x\ (125) 

We will briefly comment on the most prominent differences of Y^l and X^l- 
First, without QCD no logarithmic terms are present in Y{x), as can be seen in 
( 11251 ). They are to lowest order cancelled between the Z-penguin- and the box con- 
tribution. As a consequence the non-logarithmic term ~ x is more important here 
than in the case of X^vl- Logarithms are generated through QCD and summed to 



next-to-leading order in the renormalization group result ( 115 ). Their effect is to 
enhance the zeroth order expression by a factor of about 2.5. Recall that the func- 
tion X(xc) on the contrary is suppressed by ~ 30% through QCD effects. Nev- 
ertheless numerically X^vl still exceeds Yatl by a factor of four. Comparing the 
magnitude of the corrections in Xjvl and Yatl one realizes that QCD is much more 
important in Yatl than in X^l. The consequence is a larger relative /x-dependence 
of Ynl- In fact one finds {YNL{lGeV) - FjvL(3Ge1/))/F7VL(mc) = 43.6% whereas 
in LLA {YLilGeV) - YL^'iGeV)) /YL{m^) = 71%. Here Yl is defined to be the 
leading log function but including the numerically important term x(/i)/2, al- 
though this is not fully consistent in this approximation. The effect of a b-quark 
threshold, which we neglect, is also larger here (~ 1%) than it was for Xnl- 
The /i-dependences are illustrated in fig. 7 where the leading log result is shown 
without (Yl) and including (Y^) the non-logarithmic term x/2. Since the residual 
//-dependence turns out to be quite pronounced even in NLLA, one may start to 
worry about the vahdity of perturbation theory. At least the scale ambiguity 
is substantially reduced when going from leading logs to NLLA. We will take 
the observed variation as an estimate for the theoretical uncertainty even in this 
case and briefly discuss the consequences in the next subsection. Fortunately the 
function Y^l is less important for Kl — > fJ'~^fJ^~ than X^l is for tx^uv. 
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^M5' [GeV] 


1.2 


1.3 


1.4 


0.15 


2.78 


3.18 


3.61 


0.25 


2.96 


3.39 


3.84 


0.35 


3.10 


3.55 


4.03 



Table 5: The function Ynl for various A;^^ and rric- 

Consequently the net impact of the /i-dependence on B{Kl will turn 

out to be not as dramatic as could be expected from the discussion above. Finally 
in table 5 we give numerical values of Y^l for various nic and Ajjg with fic = 



7.3 B{Ki jj,~^fi )sD cind the Parameter g 

Expressing CKM elements through Wolfenstein parameters we find 

B{Kl ^ = 1.71 ■ 10-' A^Y\xt)[go - g]^ (126) 

where 

_ 417 F;vL 

Formula ( |126| ) together with Y{xt) and Yatl in (|112|) and ( p.l5| ) respectively gives 
B{Kl —>■ ^^^~)sD in the standard model with QCD effects in next-to-leading 
order. It is the generalization of the QCD calculations in |Q where only leading 
logs and some non- leading contributions have been taken into account. The 
experimental value of B{Kl — > ^i^^i~)sD determines the value of q given by 



P=-p-f f = j B{KL^f^^^i-)sn 1 

^ ° " V 1.71-10-9 A^Yixt) ^ ^ 

Similarly to in the case of n^uu, the value of fo is fully determined by 

the top contribution which has only a very weak renormalization scale ambiguity 
after the inclusion of 0{as) corrections. The main scale ambiguity is present in 
Qo whose departure from unity measures the relative importance of the charm 
contribution. In order to illustrate the scale ambiguities in the context of Kl — >■ 
/i"*"/!" we show in table 6 fo, Qo, q and B{Kl — > fi'^fJ'~)sD and their uncertainties. 
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NLLA (LLA) 


ro 


Qo 


Q 


B{Kl ^ /i+/i-)5D/10-9 


/^c(t) = rnc(t) 


1.12 (1.16) 


1.21 (1.20) 


0.09 (0.04) 


1.00 (0.91) 


mill 


1.12 (1.10) 


1.14 (1.10) 


0.02 (-0.14) 


0.88 (0.67) 


max 


1.12 (1.26) 


1.24 (1.27) 


0.12 (0.14) 


1.05 (1.09) 




0.6% (13.8%) 


7.6% (13.8%) 




16.8% (45.8%) 



Table 6: Theoretical uncertainties in the phenomenology of Kl — > /i+yU in LLA 
and NLLA. 

Here the same input is used as for table 4 (sect. 6.3.). To obtain the numbers in 
the last column we have set g = 0.09, corresponding to B{Kl 1^'^ IJ'~)sd = 10~^ 
in NLLA. 

In the case of q the uncertainty Af? = 0.28 present in LLA is reduced to /\q = 0.10 
in NLLA. Similarly the 46% uncertainty in B{Kl — > fi^fi~)sD is reduced to 17%. 
This theoretical uncertainty in the short- distance part of Kl f^^fJ'~ left over 
even after inclusion of next-to-leading order corrections together with the poorly 
known long distance contributions makes this decay unfortunately less suitable 
for the determination of standard model parameters than tt^uu. 

8 Summary 

In this paper we have given a detailed analysis of the theoretical prediction for 
the rare decay n^uu in the standard model. We have performed the first 

calculation of QCD effects in this process at the next-to-leading order. This re- 
quired two different approaches, the inclusion of the full 0{as) correction to all 
orders in mt/M]^ for the top-quark contribution [|^, |15[ and a two-loop renor- 
malization group calculation in the charm quark case performed here. Compared 
to the previous leading-log estimates ^, Q the main benefits of our analysis 
are: 

• The scale ambiguity in the top quark mass can be essentially removed. 

• The scale dependence is substantially reduced in the charm sector. Conse- 
quently the uncertainty due to the choice of mc{fi) stressed in [^, || is also 
reduced automatically. 
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• The next-to-leading renormalization group calculation allows a meaningful 
use of for the charm contribution. 

• The non-logarithmic mass terms ~ Xc and the effect of the tau-lepton mass 
can be incorporated consistently. 

• The renormalization scheme dependence in QCD corrections to non-loga- 
rithmic terms addressed in |^ can be avoided. 

The consequence of these improvements is a considerably reduced uncertainty of 
the theoretical prediction for B{K^ — > tt+z/z/) and in the related determination 
of CKM parameters from this decay. 

As a byproduct of the analysis of n^uu the short-distance part of Kl 

fi^fj^~ has been obtained in NLLA. In this case the scale uncertainties are also 
reduced considerably compared to the LLA, but the residual ambiguity is larger 
than in tt^uu. Still B{Kl f^^fJ'~)sD could be used to constrain the 

Wolfenstein parameter g with improved accuracy in NLLA. Unfortunately how- 
ever the presence of the poorly known long distance contribution to Kl — ^ ^'^ ^~ 
makes this decay less suitable for a determination of standard model parameters 
than tt^uu. 

71^ uu could be used in particular to constrain and eventually determine 
Vtd- B° — B° mixing, also sensitive to Vu, is probably less suited for this purpose 
because it involves a hadronic matrix element, which is hard to calculate reliably. 
By contrast, no further theoretical uncertainty is present in n^uu besides 

the residual uncertainty due to the truncation of the (RG improved) perturbation 
series, which illustrates the clean nature of this decay. For the determination of 
I Vtd I this ambiguity amounts to roughly 7%, compared to about 30% in LLA, for 
typical parameter values. The remaining uncertainty is almost entirely due to 
the charm sector. 

We have also analyzed the impact of next-to-leading corrections to n^uu 
on the determination of the parameters {g, rj), relevant for the unitarity triangle. 
The inclusion of these corrections considerably reduces the uncertainties in {g, rj) 
present in LLA. 

We remark in this context that, since the charm sector does not matter for the 
CP violating mode Kl 7t°i>ij, the standard model prediction in this case is ex- 
tremly accurate after the inclusion of QCD corrections. In fact this decay is then 
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essentially free of theoretical uncertainties, which makes the challenging search 
for Kl ii°vv particularly desirable. Similar comments apply to -B — Xsvv. 
For details see [p^ . 

The residual theoretical ambiguity in -n^vv is much smaller than uncer- 

tainties coming from the various input parameters involved. Clearly, to be able 
to take advantage of the clean theoretical prediction, further progress has to be 
made on this side. The most important tasks for the future are the determination 
of the top quark mass, a more accurate determination of |V^b| and above all - 
hopefully - a measurement of n^vv. 
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Figure Captions 



Fig. 1 The leading one-loop Z-penguin- and box diagrams inducing sd — > uu, 
responsible for K'^ tt'^vP. 

Fig. 2 The Z-penguin in the effective theory relevant for the charm contribution 
(W-, Z-boson integrated out). The square blob denotes operator insertion. 

Fig. 3 Gluonic corrections to the diagrams of fig. 2. Symmetric diagrams and coun- 
terterm contributions are not shown explicitly. The last diagram vanishes 
for massless, on-shell external quarks in dimensional regularization. 

Fig. 4 The box diagram in the effective theory. 

Fig. 5 Gluonic corrections to the diagram of fig. 4. 

Fig. 6 The function X representing the charm contribution to the sd — > Pu am- 
plitude (for nil — 0) and its dependence on the renormalization scale /ic, 
reflecting the theoretical uncertainty. The cases shown are: Xl (LLA), Xnl 
(NLLA) and Xq (without QCD). 

Fig. 7 The scale dependence of the function Y, relevant for sd fi/i. We show 
the results in LLA, with (Yj,) and without (Y^) non-logarithmic mass term, 
in NLLA (Ynl) and without QCD. 

Fig. 8 \Vtd\ determined from B{K^ tt'^uu) as a function of rrit = mt{mt) for 
fixed values of the necessary input parameters (see sect. 6). The narrow 
(broad) band results from a variation of the renormalization scales Hc, l-^t in 
the ranges IGeV < fic ^ 3GeV and |mt < fit ^ Sm^ in NLLA (LLA) and 
indicates the theoretical uncertainty involved in this analysis. 

Fig. 9 Theoretical uncertainties in the determination of the unitarity triangle in 
the {g,r]) plane from B{K'^ — > tt'^vD). For flxed parameter values the vertex 
of the unitarity triangle is forced to lie on a circle of radius around the 
origin. The variation of the scales fici l^'t within IGeV < /Xc < SGeV, 
lOOGeV < fit < 3006*61^ then yields the indicated ranges in LLA (full) and 
NLLA (reduced). We show the cases = 0.31, 0.45, 0.59. 

Fig. 10 The scale ambiguity in B{K'^ — > n'^uP), in analogy to flg.8. 
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